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We comment on two issues in quantum cosmology, in the context of the Wheeler-De 
Witt equation and wave function of the Universe: (i) arrow of time and interpretation 
of the wave function in the classically allowed regions; (ii) stability of an approximation 
of the Born-Oppenheimer type in classically forbidden regions of the scale factor. 
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Most of the evolution of the Universe is likely to have proceeded classically, in the sense 
that the dominant phenomenon is the classical expansion while quantum fluctuations are 
small and can be treated as perturbations. Prior to this stage, however, genuinely quantum 
phenomena almost certainly took place. Although it is not clear whether they have left any 
observable footprints, it is of interest to try to understand them, as this may shed light on 
such issues as initial conditions for classical cosmology (are inflationary initial data natural? 
how long was the inflationary epoch? is open Universe consistent with inflation?), properties 
of space-time near the cosmological singularity, origin of coupling constants, etc. 
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To describe the Universe at its quantum phase, one ultimately has to deal with full 
quantum gravity theory, well beyond the Einstein gravity. It is, however, legitimate to 
take more modest attitude and consider quantum phenomena below the Planck (or string) 
energy scale. Then the quantized Einstein gravity (plus quantized matter fields) provides 
an effective "low energy" description which must be tractable, at least in principle, within 
quantum field theory framework. Processes that should be possible to consider in this way 
are not necessarily perturbative, as the example of tunneling in quantum mechanics and 
field theory shows. 

Surprisingly, not so many phenomena are well understood even within this modest 
approach. Perhaps the most clear process is the decay of a metastable vacuum |jl|. It has 
been clarified recently || that at least in the range of parameters where the treatment 
of space-time in terms of background de Sitter metrics is reliable, the Coleman-De Luccia 
instanton indeed describes the false vacuum decay, provided the quantum fluctuations above 
the classical false vacuum are in de Sitter-invariant (conformal) vacuum state. This is in 
full accord with the results of Refs. [l], [| . It is likely that this conclusion holds also when 
the quantum properties of metrics are taken into account. Furthermore, the Hawking-Moss 
instanton || can be interpreted [T2J as a limiting case of constrained instantons that describe 
the false vacuum decay in an appropriate region of parameter space, again in agreement 
with previous analyses || ||. Hence, there emerges a coherent picture of the false vacuum 
decay with gravity effects included. 

One may try to apply laws of quantum mechanics to the Universe as a whole, and 
consider the wave function of the Universe. Although research in this direction began 
more than 30 years ago J7|, the situation here is still intriguing and controversial. The 
main purpose of this contribution is to make a few comments on this subject. Namely, we 
will discuss which analogies to ordinary quantum mechanics are likely to work in quantum 
cosmology, and which are rather misleading. We begin with quantum mechanics, and only 
then turn to the wave function of the Universe. 



2 Wave function in quantum mechanics 

To set the stage, let us consider a quantum mechanical system with two dynamical coordi- 
nates, x and y. Let the Hamiltonian be 



where 



H = Hq + Hy 



Ho = -$\ + V (x) 



Hy = ^P 2 y + ^\x)y 2 + -\{x)y A + ... 

Let us assume that the potential Vq{x) is such that the motion along the coordinate x is 
semiclassical, while the dynamics along the coordinate y can be treated in perturbation 
theory about the semiclassical motion along x. This approach is close in spirit to the 
Born-Oppenheimer approximation. We will consider solutions to the stationary Schrodinger 
equation with fixed energy E. 

Let us first discuss the dynamics in the classically allowed region of x, where E > Vq(x). 
In this region, there are two sets of solutions with the semiclassical parts of the wave 
functions equal to 

^ oc e +iS{x) (1) 

and 

^ oc e~ iS{x) 

where 

S(x) = f X alx' ^2(E - V (x')) 

These two sets of solutions correspond to motion right and left, respectively. Note that this 
interpretation is based on the fact that there exists extrinsic time t inherent in the problem: 
the complete, time-dependent wave functions are exp(— iEt + iS(x)) and exp(— iEt — iS(x)); 
the wave packets constructed out of the wave functions of these two types indeed move right 
and left, respectively, as t increases. 



Let us now consider the dynamics along the coordinate y, still using the time-independent 
Schrodinger equation in the allowed region of x. This is done for, say, right-moving system 
by writing, instead of eq. ([!]), 

#(x,y) = V(x,y)e iS{x) 

\fPx(x) 

where p x = dS/dx. To the first order in h one obtains that the time-independent Schrodinger 
equation reduces to 

dx dx y 
This can be cast into the form of time- dependent Schrodinger equation by changing variables 
from x to r related by x = x c (t), where x c (t) is the solution of the classical equation of 
motion for x in "time" r, which has energy E and obeys 

!<—■>-£ « 3 > 

After this change of variables, \l/ becomes a function of y and r and obeys the following 
equation, 

t ^lll = H y (y,p y] rMy;r) (4) 

where the explicit dependence of H y on r comes from x c (t). We see that there have emerged 
intrinsic time r which parameterizes the classical trajectory x c (t) and also the ^/-dependent 
part of the wave function. We note again that in quantum mechanics, the arrow of intrinsic 
time, which is set by the sign convention in eq.(f|), is determined by the arrow of extrinsic 
time t. 

Note also that one is free to choose any representation for operators y and p y and write, 
instead of eq. (||) , 

i^ = H y (T)\*} (5) 

To solve this equation, one may find convenient to switch to the Heisenberg representation, 
as usual. 



We now turn to the discussion of the region of x where the classical motion is forbidden 
and the system has to tunnel. To simplify formulas, we set E = in what follows. If the 
system tunnels from left to right, the dominant semiclassical wave function is 

\1> oc e~ s(x) 



where S(x) = J x dx' \/2Vq(x') and obeys the following equation, 

+ V (x) = 



1 fdS^ 2 



2 \dx / 

This equation may be formally considered as the classical Hamilton-Jacobi equation in 

Euclidean ("imaginary") time. The zero energy classical trajectory x c (r) in Euclidean time 

r obeys 

d 2 x c dV . 

~d^ = + ~dx~ {x = Xc) 

and hence 

dx c OS , . 

t~ = ir^ x = Xc > 

dr ox 
Then S(x) can be calculated as the value of the Euclidean action along this trajectory. 

To find the equation governing the dynamics along y-direction in the classically forbidden 
region of x, we again write 

W(x,y) = ^(x,y)e~ s{x) 

\fWx 

and obtain, changing variables from x to r, x = x c (t), that ^ obeys the time-dependent 
Schrodinger equation, now in Euclidean time, 

^^ = -H y (y,p y ;rMy,r) (6) 

The minus sign on the right hand side of this equation is crucial for the stability of the 
approximation we use. Indeed, the system described by eq.@ tends to de-excite, rather 
than excite, as "time" r increases, so that the part \& of the wave function remains always 
subdominant as compared to the leading semiclassical exponential. The physics behind 
this property is quite clear: we consider tunneling at fixed energy, so the de-excitation of 



fluctuations along y means the transfer of energy to the tunneling subsystem, which makes 
tunneling (exponentially) more probable. Inversely, if fluctuations along y get excited, the 
kinetic energy along x decreases, and tunneling gets suppressed stronger. 

3 Wave function of the Universe 

To discuss specific aspects of quantum cosmology, let us consider the closed Friedmann- 
Robertson-Walker Universe with the scale factor a. Let us introduce the cosmological 
constant A, minimal scalar field <f>{x) with a scalar potential V{4>) and also massless con- 
formal scalar field. We are going to treat the dynamics of the scale factor in a semiclassical 
manner; in this respect a is analogous to the variable x of the previous section. The minimal 
scalar field (as well as gravitons) will be considered within perturbation theory, so each of 
the modes 0k will be analogous to the variable y of the previous section. 

The basic equation in quantum cosmology is the Wheeler-De Witt equation, which in 
our case reads 

-\p 2 a -\a 2 + ka* + H^m = -^ (7) 

where we have set "iMpJlQit = 1 and ignored the operator ordering problems which are 
irrelevant for our discussion. Here 

d 3 x 



Ha, 



2tt 2 



^ + yW) 2 + a 4 V(0) 



is the term due to the minimal scalar field; at the classical level H^ is the energy of matter 
defined with respect to conformal time. The non-negative constant e on the right hand side 
of eq. (0) is the contribution of the conformal scalar field; the only purpose of introducing 
the latter field is to allow for non-zero e. We do not consider gravitons in what follows, as 
they are similar to the quanta of the minimal scalar field 0. 

In the spirit of the Born-Oppenheimer approximation, let us first neglect the conformal 
energy of the field 0, i.e., omit the term H^ in eq.(^). Then the Wheeler-De Witt equation 



takes the form of the time-independent Schrodinger equation in quantum mechanics of one 
generalized coordinate a with energy e and potential 

JJ(a) = -a 2 - Aa 4 

At 16A 2 e < 1, there are two classically allowed regions: at small a (0 < a 2 < [1 — 



VI - 16A 2 e]/4A) and at large a (oo > a 2 > [1 + Vl - 16A 2 e]/4A). At the classical level, the 
former region corresponds to an expanding and recollapsing Friedmann-like closed Universe, 
while the latter corresponds to the de Sitter-like behavior. As e — ► 0, the first classically 
allowed region disappears, while the second becomes exactly de Sitter. 

In between these two regions, classical evolution is impossible (if one neglects H^), and 
one has to consider classically forbidden "motion". Let us discuss classically allowed and 
classically forbidden regions separately. 

3.1 Classically allowed region: issue of arrow of time 

To be specific, let us consider classically allowed de Sitter-like region where the scale factor 
a is large. In the leading order, there are again two types of semiclassical wave functions, 

tf oc e~ lS{a) (8) 

and 

* oc e +l5(a) (9) 

where 

S(a) = r da' ^/2(e-U(a')) 

Classically, the momentum is related to the derivative of the conformal factor with respect 

to conformal time, 

da 

dr\ 



For the two semiclassical wave functions one has 

p°* = ( T f ) * 

where upper and lower signs refer to eq.(|) and eq.(H), respectively. Hence, one is tempted to 
interpret the wave functions @ and (|9]) as describing expanding and contracting Universes, 
respectively. Indeed, the Hartle-Hawking wave function || that in the allowed region is a 
superposition 

* M oce- !S(a) + e +lS(fl) (10) 

is often interpreted as describing a collapsing and re-expanding de Sitter-like Universe. 
Similar interpretation is often given to the Linde wave function ||. On the other hand, the 



tunneling wave functions Jl0|, [11], [12] which contain one wave in the allowed region, 

are often assumed to be the only ones that correspond to an expanding, but not contracting, 
Universe; this is, at least partially, the basis for the tunneling interpretation. 

An important difference with conventional quantum mechanics is, however, the absence 
of extrinsic time in quantum cosmology. Hence, the arrow of intrinsic time has yet to be 
determined. In other words, there is no a priori reason to interpret the wave functions (|8|) 
and (§) as describing expanding and contracting Universes, respectively. The sign of the 
semiclassical exponent does not by itself determine the arrow of time. 

Were the scale factor the only dynamical variable, it would be impossible to decide 
whether, say, the wave function @ corresponds to expanding or contracting Universe. If 
the matter fields (and/or gravitons) are included, this should be possible. Before discussing 



this point, let us derive the equation for the wave function describing matter [|13|, [14|, [L5 
again in the spirit of the Born-Oppenheimer approximation. 

Let us extend the wave functions fl8|) and @ to contain the dependence on the matter 
variables, 

|tf( a )\ = _Le TiS(a) ||r(a)) (11) 
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where at given a both |\l/(a)) and |^(a)) belong to the Hilbert space in which 0(x) and 
p<p("x) act. As an example, one may (but does not have to) choose the generalized coordinate 
representation; then |^(a)) becomes a function ^({0k}; a) of the Fourier components of 0. 
In the first order in h one obtains from eq.(^) 



±iy/e-U(a)^^- = H^(a)) (12) 



in complete analogy to eq. 



The arrow of time is determined now by where (at what a) and which initial conditions 
are imposed on |\l/(a)). As an example, let us assume that the initial conditions for the 
evolution in real intrinsic time are imposed at small a (at the turning point a 2 = [1 + 



\/i — 16A 2 e]/4A), and that at that point |\I/) describes smooth distribution of the scalar 
field. This type of initial data are characteristic, in particular, to the Hartle-Hawking no- 
boundary wave function. As a increases, the system will become more and more disordered, 
independently of the sign in eq. (fP2"D . With thermodynamical arrow of time, both wave 
functions (|Tl~D will describe expanding Universe. 



If, with these initial conditions, one changes variables from a to rj using 

da 



then rj increases with a, so that rj is the conformal intrinsic time, independently of the choice 
of sign in eq . (|TT|) . In the case of positive sign, eq. (|T2|) becomes the conventional Schrodinger 



equation for quantized matter in the expanding Universe, 

d\ty) 

^ = HMm (is) 

where the matter Hamiltonian depends on 77 through 0(77). On the other hand, in the case 
of negative sign one obtains "wrong sign" Schrodinger equation, 

dm 

This little problem is easily cured by considering, instead of |^), its T-conjugate, |^^}; 
if the generalized coordinate representation is chosen for |\l/), then T-conjugation is merely 



complex conjugation, \]/( T )(0k; ??) = ^*(0k;^)- The T-conjugate wave function obeys con- 
ventional Schrodinger equation, but with CP-transformed Hamiltonian. Hence, the inter- 
pretation of both wave functions (JTTJ) as describing the expanding Universe is self-consistent; 
the only peculiarity is that the wave function e +l5 |\&) corresponds to the Universe in which 
matter is CP-conjugate. In particular, we argue that both components of the Hartle- 
Hawking wave function fllHD correspond to expanding Universes. 

In more generic cases (in particular, when the matter degrees of freedom cannot be 
treated perturbatively, see, e.g., Refs. |16|, [L7|] and references therein), the situation may be 
much more complicated. Still, the arrow of time is generally expected to be one of the key 
issues in the interpretation of the wave function of the Universe. 

3.2 Classically forbidden region: issue of stability of the Born- 
Oppenheimer approximation 

We now consider the region of the scale factor that is classically forbidden in the absense of 

Hf/,, i.e., oi < a < 02, where 

2 1 =F VI - 16A 2 e 

«1,2 - ^ 

If Hf/) is switched off, there are two semiclassical solutions to the Wheeler-De Witt equation, 

* oc e" 5(a) (14) 

and 

^oce +s(a) (15) 

where 

S(a) = f a da' y/2(U(a') - e) 

is defined in such a way that it always increases at large a. The wave function flTjj) decays 
as a increases, so it may be interpreted as describing tunneling from classically allowed 
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Friedmann region to de Sitter-like one. It is convenient to introduce Euclidean conformal 
time parameter r and consider Euclidean trajectory a c (r) obeying 

da c dS , . 

At e = the Euclidean four-geometry corresponding to this trajectory is a four-sphere, the 
standard de Sitter instanton. 

Let us now turn on the scalar field Hamiltonian H,p, and try to apply the procedure of 



the Born-Oppenheimer type. We write, instead of eq.fll4|), for the wave function decaying 
at large a, 

Ma)) = ^e- 5(a) |^(a)) 



and obtain in the first order in h that |^(a)) obeys the "wrong sign" Euclidean Schrodinger 

0|%)) 



equation |12 



+H+(t)\9(t)) (16) 

where the change of variables from a to r with a = a c (r) has been performed. The sign on 
the right hand side of eq-flTB]) is opposite to that appearing in the usual quantum mechanics, 
eq.(0), and is directly related to the sign of p^-term in the Wheeler-De Witt equation (|?p. 

The "wrong" sign in eq.flTop implies that the approximation we use is in fact unstable, 
if generic "initial" conditions are imposed at small a, say, at a — a\. Note that imposing 
initial conditions in this way is natural if one interprets the wave function decaying at large 
a as describing tunneling from small to large a. The formal reason for the instability of the 
approximation is that the degrees of freedom of the scalar field get excited as a increases in 
the forbidden region. The rate at which this excitation occurs is generically high [1^], and 



the approximation breaks down well before a gets close to the second turning point 02- 

In the path integral framework, breaking of the Born-Oppenheimer-type approxima- 
tion for the wave function decaying at large a is also manifest fll8| . This wave function 
corresponds to the Euclidean path integral with "wrong" sign of the action, 

f Dg D(j) e +5[ ^ ] 
11 



The instanton action then gives the factor e Sinst , but the integral over (and gravitons) 
diverges. 

The physics behind this instability is that tunneling of a Universe filled with matter is 
exponentially more probable as compared to empty Universe. Hence, the matter degrees 
of freedom tend to get excited in the forbidden region, thus making tunneling easier. Note 
that this property is peculiar to quantum cosmology: in quantum mechanics the situation 
is opposite, as we discussed in the previous section. 

There are exceptional cases in which matter degrees of freedom do not get excited in the 
forbidden region, e.g., because of symmetry. In our model this would be the case if e = 
and the scalar field was in the de Sitter-invariant state, cf. Ref. |1|. Such cases do not 
seem generic, however. 

Breaking of the Born-Oppenheimer approximation does not necessarily mean that 
tunneling-like transitions from small a to large a with generic state of matter at small 
a do not make sense. Rather, it is the semiclassical expansion that does not work in this 
case, so the state of the Universe after the transition may be quite unusual. Presently, 
neither the properties of this state, nor the properties of the wave function in the forbidden 
region are understood (except for special cases mentioned above). 

The situation is different for the wave functions increasing towards large a, eq . fllTf) . 
In that case the matter wave function obeys the usual Euclidean Schrodinger equation, 
<9|\I/(t))/<9t = —H^(t)\^(t)), where r is still assumed to increase with a. Hence, it is 
possible to impose fairly general initial conditions at small a, and the approximation will not 
break down. In particular, the Hartle-Hawking wave function is a legitimate approximate 
solution to the Wheeler-De Witt equation in the forbidden region. This is in accord with 
the path-integral treatment: the increasing wave function ( |I5|) corresponds to the standard 
sign of the Euclidean action in the path integral. 

The non-semiclassical behavior of the tunneling wave functions, signalled by the instabil- 
ity of the Born-Oppenheimer-type approximation, is a special, and potentially interesting, 
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feature of quantum cosmology. It is a challenging technical problem to develop techniques 
adequate to this situation. It is not excluded also that the properties of the tunneling 
wave functions are rich and complex, and that understanding them may shed light on the 
beginning of our Universe. 

The author is indebted to A. Albrecht, J. Goldstone, N. Turok, W. Unruh and A. 
Vilenkin for helpful discussions. 
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